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The 'Powerful Pocketful': an Electronic 
Calculator Challenges the Slide Rule 

This nine-ounce battery-powered scientific calcu- 
lator small enough to fit in a shirt pocket, has 
logarithmic, trigonometric, and exponential functions 
and computes answers to 10 significant digits 

By Thomas M. Whitney, France Rode, and Chung C. Tung 



WHEN AN ENGINEER OR SCIENTIST 
NEEDS A QUICK ANSWER to a problem 
that requires multiplication, division, or transcen- 
dental functions, he usually reaches for his ever- 
present slide rule. Before long, 
however, that faithful 'slip stick' 
may find itself retired. There's 
now an electronic pocket cal- 
culator that produces those an- 
swers more easily, more quickly, 
and much more accurately. 

Despite its small size, the new 
HP-35 is a powerful scientific 
calculator. The initial goals set 
for its design were to build a 
shirt-pocket-sized scientific cal- 
culator with four-hour operation 
from rechargeable batteries at a 
cost any laboratory and many 
individuals could easily justify. 
The resulting nine-ounce prod- 
uct surprises even many who are 
acquainted with what today's 
large-scale integrated circuits 
can achieve. 

The HP-35 has basically the 
same functions and accuracy of 
other HP calculators, and it is 
portable. It is a close cousin to 
the many four-function elec- 
tronic calculators which have 
appeared in recent years, ini- 
tially from Japan and now from 
many U.S. manufacturers. However, three features 
set the HP-35 apart from four-function calculators. 
First, none of the four-function calculators has 
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transcendental functions (that is, trigonometric, log- 
arithmic, exponential) or even square root. Second, 
the HP-35 has a full two-hundred-decade range, al- 
lowing numbers from 10"" to 9.999999999 X lO* 0 " 
to be represented in scientific 
notation. Third, the HP-35 has 
five registers for storing con- 
stants and results instead of just 
one or two, and four of these 
registers are arranged to form 
an operational stack, a feature 
found in some computers but 
rarely in calculators (see box, 
page 5). On page 7 are a few 
examples of the complex prob- 
lems that can be solved with the 
HP-35. 




Data Entry 

The photograph of the calcu- 
lator on this page shows how 
the keys are arranged. Numbers 
enter the display, which is also 
called the X register, from left 
to right exactly as the keys are 
pressed. Entry is entirely free- 
field, that is, digits will be dis- 
played exactly as they are en- 
tered, including leading or trail- 
ing zeros. 

The enter-exponent key, EEX. 
is used for entering numbers in 
scientific notation. For example, 
the number 612,000 may be entered as 6.12 EEX 5 or 
612 EEX 3 or .0612 EEX 7 as well as 612000. The 
change-sign key, CHS, changes the sign of the man- 
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tissa or. if pressed immediately following EEX, the 
sign of the exponent. Mistakes during data entry 
can be corrected by use of the clear x key. CLx. A 
special key is provided for entering jr. 

Display 

The display consists of 15 seven-segment-plus- 
decimal-point light-emitting-diode (LED) numerals. 
Answers between 10'" and 10" : will always be dis- 
played as floating-point numbers with the decimal 
point properly located and the exponent field blank. 
Outside this range the HP-35 displays the answer in 
scientific notation with the decimal point to the 
right of the first significant digit and the proper 
power of 10 showing at the far right of the display. 
To make the display more readable, a separate digit 
position is provided for the decimal point. 

The display is always left justified with trailing 
zeros suppressed. For instance, the answer to 3 -f- 4 
appears as .75 instead of 0.750000000. 



Answers greater than 10" 10 (overflow) are dis- 
played as 9.999999999 X 10", and answers smaller 
than 10"° are returned as zeros. These are the 
HP-35's 'closest' answers. Improper operations, 
such as division by zero or the square root of a 
negative number, cause a blinking zero to appear. 

Single-Operand Function Keys 

The function keys operate on the number dis- 
played in the X register, replacing it with the func- 
tion of that argument. The trigonometric functions, 
sin. cos. tan, arc sin. arc cos. and arc tan operate in 
degrees only. For inverse trigonometric functions 
the function key is preceded by the arc key. The 
angle obtained will be the principal value. 

The functions log x and In x compute the log- 
arithm of x to base 10 and base e respectively. The 
exponential function e* is also provided. Other sin- 
gle-operand functions are the square root of x. \ r x, 
and the reciprocal of x, 1/x. 

Arithmetic Keys: Two-Operand Functions 

The arithmetic functions. +, — , X, and 4-, and 
the power function x T . operate on the X and Y regis- 
ters, with the answer appearing in X. Numbers are 
copied from X into Y by use of the ENTERt key. 
Thus to raise 2 to the 7th power the key sequence 
is 7, ENTERt, 2, x r and the answer 128 is displayed 
immediately after x r is pressed. This is a general 
principle: when a key is pressed, the correspond- 
ing operation is performed immediately. 

Although raising x to a power can be accom- 
plished via the formula x T = e-" * and in fact is 
done this way internally, the single-keystroke op- 
eration is much more convenient. The power, y. can 
be any positive or negative integer or fraction, while 
x must be positive. 

Registers and Control Keys 

Five registers are available to the user. Four of 
these, called X. Y. Z. and T, form an operational 
stack. The fifth, S, is for constant storage. (For 
clarity, capital letters refer to the registers and 
lower-case letters refer to the contents of the 
registers.) 

Five keys are used to transfer data between regis- 
ters. The ENTERt key pushes the stack up. that is. 
x is copied into Y, y into Z. z into T, and t is lost. 
This key is used as a separator between consecu- 
tive data entries, such as in the power example 
just described. Roll down. Rl. is used to view the 
contents of X, Y. Z. and T. Four consecutive opera- 
tions of this key return the registers to the original 
state with no loss of data. 
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"> 3 in display (X register) 

' 3 duplicated into Y register by ENTER? 

1 4 in display. 

Product (12) appears in X and stack drops. 
1 Automatic ENTERt pushes 12 into Y. Display shows 5. 
I ENTER? pushes y into I, x into Y. > is unchanged. 

6 in display. 
1 Product (30) appears in X and stack drops. 
' Sum (42) appears in X and stack drops. 



Fig. 1. HP-35 Pocket Calculator has a tour-register opera- 
tional stack (last-in-tirst-out memory). Here's how the 
stack works to solve (3 X A) + (5 X 6). Answers appear in 
display register X. in lloating-point or scientitic notation, 
to 10 significant digits. 

The contents of X and Y can be exchanged using 
jTy; this is useful if the operands have been entered 
in the improper order for x T , — , or -K The constant 
storage location, S. is accessed via the two keys 
store, STO. and recall, RCL. 

The operational stack is automatically pushed up 
for data entries following any operation other than 
ENTERt . STO. and CLx. This saves many uses of the 
ENTERt key. The stack automatically drops down 
following any two-operand function ( + , — , X, 
x r ). For example, to do the problem (3X4) + (5X6) 
the key sequence is: 3, ENTERt. 4, X. 5, ENTERt. 
6, X, + (see Fig. 1). ENTERt is unnecessary be- 
tween X and 5 because the answer to the first term, 
12, is automatically transferred to Y when the 5 is 
entered. Also, no RVis necessary after the second X 
since the first term, 12. is automatically transferred 
from Z to Y after the multiplication. 

The clear-all key. CLR, clears all registers includ- 
ing storage. Initial turn-on of the calculator has the 
same effect. 

On the back of the calculator is an instruction 
panel that provides the user with quick answers to 
the most commonly asked 'how-to-do-it* questions 
(Fig. 2). The panel also shows an example of a prob- 
lem solution. 

System Organization 

Now let's go inside and see what makes it work. 
The HP-35 contains five MOS/LSI (metal-oxide- 



Fig. 2. Instruction panel on back ot calculator answers 
most Irequently asked questions. 

semiconductor/large-scale-integration) circuits: 
three read-only-memories (ROMs), an arithmetic 
and register circuit (A&R), and a control and timing 
circuit (C&T). The logic design was done by HP and 
the circuits were developed and manufactured 
by two outside vendors. Three custom bipolar cir- 
cuits are manufactured by HP's Santa Clara Divi- 
sion: a two-phase clock driver, an LED anode 
driver/clock generator, and an LED cathode driver. 
Fig. 3 is a block diagram of the calculator. 

The HP-35 is assembled on two printed circuit 
boards (see Fig. 4). The upper board contains the 
display and drivers and the keyboard. The lower 
and smaller board has all the MOS logic, the clock 
driver, and the power supply. 

The calculator is organized on a digit-serial, bit- 
serial basis. This organization minimizes the num- 
ber of connections on each chip and between chips, 
thereby saving area and cost and improving reliabil- 
ity. Each word consists of 14 binary-coded-decimal 
digits, so each word is 56 bits long. Ten of the 14 
digits are allocated to the mantissa, one to the 
mantissa sign, two to the exponent, and one to the 
exponent sign. 

Three main bus lines connect the MOS circuits. 
One carries a word synchronization signal (SYNC) 
generated by a 56-state counter on the control and 
timing chip. On another bus. instructions (Is) are 
transmitted serially from the ROMs to the control 
and timing chip or to the arithmetic and register 
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Operational Stacks and 
Reverse Polish Notation 

In 1951, Jan Lukasiewicz' book on formal logic first demon- 
strated that arbitrary expressions could be specified unam- 
biguously without parentheses by placing operators im- 
mediately before or after their operands. For example, the 
expression 

fa + b) x (c - d) 

is specified in operator prefix notation as 

X + ab - cd 

which may be read as multiply the sum a plus b by the dif- 
ference c minus d. Similarly, the expression can be speci- 
fied in operator postfix notation as 

ab + cd - X 

with the same meaning. In honor of Lukasiewicz, prefix and 
postfix notation became widely known as Polish and reverse 
Polish, respectively. 

During the following decade the merits of reverse Polish 
notation were studied and two simplifications in the execu- 
tion of computer arithmetic were discovered. First, as re- 
verse Polish notation is scanned from left to right, every 
operator that is encountered may be executed immediately 
This is In contrast to notation with parentheses where the 
execution of operators must be delayed In the above ex- 
ample, (a + b) x (c — d). the multiply must wait until 
(c — d) is evaluated This requires additional memory and 
bookkeeping Second, if a stack (that is, a last-in-first-out 
memory) is used to store operands as a reverse Polish ex- 
pression is evaluated, the operands that an operator re- 
quires are always at the bottom of the stack (last operands 
entered) For (a + b) X (c - d), the reverse Polish, 
ab + cd — x . is evaluated as follows. 
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Reverse 
Polish 


Enter a 


Enter b 
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Enter c 


Enter a 




X 



These properties have made the notation a valuable tool 
in the computer industry. All modern computer compilers for 
languages such as FORTRAN and ALGOL convert state- 
ments to reverse Polish in some form before producing a 
program that can be executed Some computer manufac- 
turers have even designed their machines with special in- 
structions to perform stack operations to facilitate execution 
of reverse Polish However, the HP-35 is the first scientific 
calculator to fully exploit the advantages of reverse Polish 
and automatic stack operations to provide user convenience 
seldom found in calculators. 



chip. The third bus signal, called word select (WS), 
is a gating signal generated on the C&T chip or by 
the ROMs: it enables the arithmetic unit for a por- 
tion of a word time, thereby allowing operations 
on only part of a number, such as the mantissa or 
the exponent. 



Control and Timing Circuit 

The control and timing (C&T) circuit performs 
the major nonarithmetic, or housekeeping functions 
in the calculator. These include interrogating the 
keyboard, keeping track of the status of the system, 
synchronizing the system, and modifying instruc- 
tion addresses. 

The keyboard is arranged as a five-column, eight- 
row matrix. It is scanned continuously by the C&T 
chip. When contact is made between a row and a 
column by pressing a key, a code corresponding to 
that row and column is transmitted over the L line 
to the read-only memory (ROM). This code is the 
starling address of a program in ROM to service that 
key. Key bounce and lockout are handled by pro- 
grammed delays. 

In all digital systems, status bits or flags are used 
to keep track of past events. In the HP-35 there are 
twelve status bits, all located on the C&T chip. They 
can be set, reset, or interrogated by microinstruc- 
tions. 

ROM addresses are updated on the C&T chip and 
sent serially to the ROMs over the It bus. During 
execution of a branch instruction, the appropriate 
signal arithmetic carry or status bit — is tested to 
determine whether the incremented address or the 
branch address should be selected next. 

A powerful feature of the serial organization is 
the ability to operate on just a single digit of a num- 
ber as it flows through the arithmetic unit. On the 
C&T chip are a pointer register and a word-select 
circuit which issue a word-select signal (WS) cor- 
responding to the time slot being operated upon. 
The value of the pointer register can be set, incre- 
mented, and decremented by microinstructions. 

Read-Only Memory 

Preprogrammed mathematical routines are stored 
in three ROM chips, each of which contains 256 
instructions of 10 bits each. A specific select code 
is assigned to each ROM chip. Only one of the three 
ROM chips is used at any time. When a ROM selec- 
tion instruction is issued, a decoder inside each 
ROM checks the select-code field of the instruction. 
In case of a match, the selected ROM turns on. Un- 
selected ROMs turn off. 

A timing circuit on each ROM is synchronized to 
the SYNC signal issued by the C&T chip as the cal- 
culator is turned on. This circuit then keeps the 
ROM chip running synchronously with the rest of 
the system. 

A ROM address register on each ROM chip re- 
ceives the address sent out by the C&T chip. The 
corresponding instruction is placed on the instruc- 
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Three custom HP-developed bipolar IC's are manufactured by HP. SYNC. Is. and WS are 
three main bus lines connecting the MOS circuits. 

(ion line. Is. provided the ROM chip is turned on. 

The ROM chip also issues word-select signals 
for certain classes of instructions. 

Arithmetic and Register Circuit 

The arithmetic and register circuit executes in- 
structions coming in bit-serially on the Is line. Most 
arithmetic instructions must be enabled by WS, the 
word-select signal. Data to be displayed is sent to 
the LED anode drivers on five lines, and one carry 
line transfers carry information back to the C&T 
chip. The BCD output is bidirectional and can carry 
digits into and out of the A&R chip. 

The A&R circuit is divided into five areas: in- 
struction storage and decoding circuits, a timing 
circuit, seven 56-bit registers, an adder-subtractor, 
and a display decoder. 

Three of the registers are working registers. One 
of these and three of the remaining four registers 
form the four-register stack. The seventh register 
is an independent register for constant storage. 
There are numerous interconnections between reg- 
isters to allow for such instructions as exchange, 
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How the HP-35 Compares with the Slide Rule 



These comparisons were made by engineers who are no! 
only highly proficient in slide rule calculation, but were also 
familiar with the operation of the HP-35. Thus, the solution 
times should not be taken as typical They do. however, 
serve to indicate the relative time advantage of the HP-35 
and to point up the still more significant advantage of its 
accuracy. 

PROBLEM U COLLECTION SOLID ANGLE FROM A POINT 
SOURCE. 
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HP-35 SOLUTION: 

2.5 Gj 10.3iZXI]La31[±JSlBSGDlLt]2CiD'3:GEI —- 

.1772825509 

SLIDE RULE SOLUTION: = 

TIME ON HP-35; 20 seconds with answer to ten significant 
digits. 

TIME ON SLIDE RULE: 3 minutes, 15 seconds with answer 
to three decimal places. 

PROBLEM 2: GREAT CIRCLE DISTANCE BETWEEN SAN 
FRANCISCO AND MIAMI. 
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HP-35 SOLUTION: 

52.4 64 3 ~ 'jl 52.4 ••• 64.3 i" _? 42.3 «- « ? •" «- 
60 ' — 2254.093016 



SLIDE RULE SOLUTION: 



tiff «M> 



TIME ON HP-35 65 seconds with answer to ten significant 
digits. 

TIME ON SLIDE RULE: 5 minutes with answer to four sig- 
nificant digits. 

PROBLEM 3. pH OF A BUFFER SOLUTION. 
HP-35 SOLUTION: 

7.2 1 CE 1 0 32 .03dD 1 E 2. 1 6 ZLi 1 0 .0087 [E ES [«j.03GD 

11JElOCaE.0087CD7^1S103SCSSSlEJ(53--- 
-7.47877778 

SLIDE RULE SOLUTION: -2*1- 

TIME ON HP-35: 65 seconds with answer to len significant 
digits. 

TIME ON SLIDE RULE: 5 minutes with answer to three sig- 
nificant digits. 



transfer, rotate stack, and so on. An advantage of 
the bit-serial structure is that interconnections re- 
quire only one gate per line. 

Transfers into or out of the stack or the constant 
register are always whole-word transfers. All other 
arithmetic instructions are controlled by the word 
select signal, WS. Thus it's possible to interchange 
only the exponent fields of two registers, or to add 
any two corresponding digits of two registers. 

The adder-subtractor computes the sum or dif- 
ference of two decimal numbers. It has two data 
inputs, storage for carry or borrow, and sum and 
carry/borrow outputs. For the first three clock 
times, the addition is strictly binary. At the fourth 
clock time the binary sum is checked, and if the 
answer is more than 1001 (nine), then the sum is 
corrected to decimal by adding 0110 (six). The re- 
sult is then entered into the last four bits of the re- 
ceiving register and the carry is stored. A similar 
correction is done for subtraction. Carry informa- 
tion is always transmitted, but is recorded by the 
control and timing chip only at the last bit time of 
the word-select signal. 



Simulation and Test of the MOS Circuits 

In designing elaborate integrated circuits like the 
C&T. A&R, and ROM chips, two questions that have 
to be answered at the very beginning are: How is 
the design to be checked? and How is the final in- 
tegrated circuit to be tested? 

The first question has two answers. One is to 
build a breadboard and compare its operation with 
the desired operation. A second answer is a com- 
puter simulation of the circuit. When the MOS cir- 
cuits (C&T, A&R. and ROMs) for the HP-35 were 
being designed, the computer simulation approach 
was chosen over a TTL or MOS breadboard. It was 
felt that the hardware breadboard wouldn't be an 
exact model of the final circuits anyway, and two 
or three months of development time could be saved 
by computer simulation because people could work 
in parallel rather than serially on a breadboard. 

A general-purpose simulation program, HP- 
DABEL. had just been developed by Jim Duley of 
HP Laboratories. This was used to check out each 
gate, each circuit, each chip, and finally all the chips 
together. Each MOS circuit is designed as a network 
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of gates and delay elements. For each gate output 
an algebraic equation was written as function of 
the inputs to the gate. This produced a large set 
of algebraic equations to be evaluated every clock 
time. A printout was available so the operation of 
any of the gates or delay outputs could be observed, 
as if with an oscilloscope probe. In this respect the 
computer simulation was much better than a hard- 
ware breadboard. 

Because of the large number of equations to be 
evaluated each clock time the general-purpose 
simulation program was too slow to use for evaluat- 
ing the algorithms implemented in the ROMs. For 
this a higher-level simulation was used, so only the 
input/output functions of each subsystem had to be 
specified. This was fast enough that all the algo- 
rithms could be checked, even the transcendental 
functions. If anything went wrong it was always 
possible to stop the program and step through it 
until the trouble was found. Correcting a problem 
was a simple matter of changing a punched card 
or two. These are advantages a hardware bread- 
board doesn't have. 

The simulation approach proved very successful. 
It saved a lot of time not only in logic design, but 
also in generating the test patterns to be used for 
testing the final integrated circuits. After a simu- 
lation is running successfully a pattern for each in- 
put is specified such that virtually every circuit 
element will be exercised. By running the program 
and recording all the inputs and outputs a complete 
test pattern is generated. This is recorded on tape, 
ready for final test of the IC. 




Fig. 5. Light-emitting-diode cluster was specially devel- 
oped lor the HP-35. Magnifying lenses are built in. 



Display and Drivers 

It was apparent early in the HP-35 planning that 
new display techniques would be required. Existing 
light-emitting-diode products used too much power 
and cost too much. HP Associates developed a mag- 
nified five-digit cluster which saves both power and 
cost and is packaged in a convenient 14-pin pack- 
age (Fig. 5). Each digit has a spherical lens molded 
in the plastic over it. A slight reduction in viewing 
angle results, but for the handheld calculator this 
is not a problem. 

LED's are more efficient if they are pulsed at a 
low duty cycle rather than driven by a dc source. 
In the HP-35, energy is stored in inductors and 
dumped into the light-emitting diodes. This drive 
technique allows a high degree of multiplexing; the 
digits are scanned one at a time, one segment at a 
time. 

Customized bipolar anode and cathode driver 
circuits incorporating the required features were 
developed and are manufactured by HP. The anode 
driver generates the two-phase system clock and 
the segment (anode) drive signals, decodes the data 
from the arithmetic and register chip and inserts 
the decimal point, sends shift signals to the other 
axis of the multiplex circuitry (the cathode driver), 
and senses low battery voltage to turn on all deci- 
mal points as a warning that about 15 minutes of 
operating time remains. The cathode driver contains 
a 15-position shift register which is incremented 
for each digit position. 

Keyboard 

Requirements for the HP-35's keyboard were 
particularly difficult. The keyboard had to be re- 
liable, inexpensive, and low-profile, and have a 
good 'feel'. The solution is based on the 'oilcan' or 
'cricket' principle, that is, curved metal restrained 
at the edges can have two stable states. The larger 
board in Fig. 4 shows the etched metal keyboard 
strips running horizontally. At each key location the 
metal is raised to form a hump over a printed-circuit 
trace running underneath. Depressing a key snaps 
the metal down to make contact with the trace. The 
keys have a definite 'fall-away' or 'over-center' feel 
so that there is no question when electrical contact 
is complete. 
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Algorithms and Accuracy 
in the HP-35 

A lot goes on in that little machine when it's com- 
puting a transcendental function. 

By David S. Cochran 



THE CHOICE OF ALGORITHMS FOR THE 
HP-35 received considerable thought. Power 
series, polynominal expansions, continued frac- 
tions, and Chebyshev polynominals were all con- 
sidered for the transcendental functions. All were 
too slow because of the number of multiplications 
and divisions required. The generalized algorithm 
that best suited the requirements of speed and pro- 
gramming efficiency for the HP-35 was an iterative 
pseudo-division and pseudo-multiplication method 
first described in 1624 by Henry Briggs in 'Arith- 
metica Logarithmica' and later by Voider' and Meg- 
gitt\ This is the same type of algorithm that was 
used in previous HP calculators. 

An estimate of program execution times was 
made, and it became apparent that, by using a bit- 
serial data word structure, circuit economies could 
be achieved without exceeding a one-second com- 
putation time for any function. Furthermore, the 
instruction address and instruction word could be 
bit-serial, too. 

The complexity of the algorithms made multilevel 
programming a necessity. This meant the calculator 
had to have subroutine capability, as well as special 
flags to indicate the status and separations of var- 
ious programs. In the HP-35, interrogation and 
branching on flag bits — or on arithmetic carry or 
borrow— are done by a separate instruction instead 
of having this capability contained as part of each 
instruction. This affords a great reduction in in- 
struction word length with only a slight decrease 
in speed. 

The arithmetic instruction set was designed spe- 
cifically for a decimal transcendental-function cal- 
culator. The basic arithmetic operations are 
performed by a 10's complement adder-subtractor 
which has data paths to three of the registers that 



are used as working storage. Partial word desig- 
nators (word select) are part of the instruction word 
to allow operating on only part of a number — for 
example, the mantissa or the exponent field. 

Sine Algorithm 

The sine routine illustrates the complexities of 
programming a sophisticated calculator. First, de- 
grees are converted to radians by multiplying by 
2:7/360. Then integer circles are removed by repeat- 
edly subtracting 2- from the absolute value of the 
argument until the result is less than 2-. If the result 
is negative, 2- is added to make it positive. Further 
prescaling to the first quadrant isn't required. The 
resulting angle is resolved by repeatedly subtracting 
tan"' 1 and counting until overdraft, then restoring, 
repeatedly subtracting tan"' 0.1 and counting until 
overdraft, and so on. This is very similar to division 
with a changing divisor. Next the resulting pseudo- 
quotient is used as a multiplier. Beginning with an 
X vector of 1 and a Y vector of 0 a fraction of X is 
added to Y and a fraction of Y is subtracted from X 
for the number of times indicated by each multiplier 
digit. The fraction is a negative power of 10 cor- 
responding to that digit position. The equations of 
the algorithm are: 

pseudo-division |fl n .i = 0 n — tan"' k 

pseudo- (Xn., = X n — Y n k 

multiplication |Y., = Yn + X„ k 

X, = 1. Y„ = 0 

k = 10 ' j = 0,1.2 . . . 

The pseudo-multiplication algorithm is similar to 
multiplication except that product and multiplicand 
are interchanged within each iteration. It is equiva- 
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Fig. 1. Accuracy ot exponential function in HP-35 Calculator Error bound is approximately 
Se'. where S Is the error due to prescaiing and the algorithm itseil S is estimated to be 
equivalent to one count in the tenth significant digit ot the argument x. 



lent to a rotation of axes. The resultant Y and X 
vectors are proportional to the sine and cosine 
respectively. The constant of proportionality arises 
because the axis rotation is by large increments 
and therefore produces a stretching of the unit 
circle. Since this constant is the same for both sine 
and cosine their ratio is identically equal to the 
tangent. The signs of each are preserved. The sine 
is derived from the tangent by the relationship 
tan 6 

Sin "- p + tan*)'"' 

Accuracy and Resolution 

Determination of the accuracy of the HP-35 is as 
complex as its algorithms. The calculator has in- 
ternal roundoff in the 11th place. In add, subtract, 
multiply, divide, and square root calculations the 
accuracy is rtVi count in the 10th digit. In calcu- 
l.'iting the transcendental functions many of these 
elementary calculations are performed with the 
roundoff error accumulating. In the sine compu- 
tation there is a divide, a multiply, and a subtract 
in the prescale operation, and there are two di- 
vides, a multiply, an addition, and a square root 
in the post-computation. Roundoff errors in these 
calculations must be added to the error of the basic 
algorithm to get the total error. 

Accuracy and resolution are sometimes in con- 
flict; for example, the subtraction of .9999999999 
from 1.0 yields only one digit of significance. This 
becomes very important, for example, in computa- 
tions of the cosines of angles very close to 90°. The 
cosine of 89.9° would be determined more accu- 
rately by finding the sine of 0.1 °. Similarly, the sine 
of 10" wastes all ten digits of significance in speci- 
fying the input angle, because all integer circles will 
be discarded. 

For many functions there is no simple exact ex- 
pression for the error. The exponential function is 
a good example. Let S be the accumulated prescai- 



ing error and computational error in the algorithm, 
referred to the input argument x. Then for 8<<1, 

e >-« - e' = e*e" - e* = e'(e* - 1) = 8e*. 

Fig. 1 shows the error bound for the exponential 
function for various arguments, assuming that 8 is 
equivalent to one count in the tenth significant 
digit of x. S 
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Packaging the Pocket Calculator 

The industrial design of the HP-35 was of primary 
importance, often taking precedence over electrical 
considerations 

By Edward T. Liljenwall 



THE INDUSTRIAL DESIGN OF THE HP-35 was 
unusual not only for Hewlett-Packard, but for 
the electronics industry in general. Usually, the 
mechanical and electrical components of a product 
are determined before the exterior is designed. The 
HP-35 took the opposite approach. 

Since the calculator was 
to be pocket-sized, size was 
the overriding constraint 
on the design. In addition 
to size three other param- 
eters were established. 
The calculator would have 
thirty to thirty-five keys, 
contain two or three bat- 
teries, and have a twelve- 
to-fifteen-digit LED dis- 
play. 

The industrial design be- 
gan with an investigation 
of keyboard, packaging, 
and overall shape concepts. 
Several basic form factors 
were studied using sketches 
and simple three-dimen- 
sional models. The models 
were particularly valuable 
at this stage of develop- 
ment. They allowed a good 
evaluation of the shapes 
and sizes being considered. 

Once the preferred direction had been established 
a detailed model was built. The model was well 
accepted and approval was given to develop the 
concept into a reality. 

Only a general idea of the electronic design was 




known at this point. Designing and packaging all 
necessary electrical and mechanical components 
into the tiny product became a tremendous chal- 
lenge lor electrical, mechanical, and industrial de- 
signers alike. 

From a human-engineering standpoint, the key- 
board was the most criti- 
cal area of the design. The 
problem was to place thir- 
ty-five keys in an area ap- 
proximately 2V2 inches by 
4V2 inches and retain the 
ability to operate the keys 
without striking more than 
one at a time. It became 
apparent that the industry 
standard of V-i-inch center- 
to-center key spacing could 
not be maintained. A suc- 
cessful compromise was to 
use n /i6-inch center-to-cen- 
ter spacing for the numeric 
keys, and Vi-inch spacing 
for all others. This was 
made possible by reduc- 
ing the size of each key. 
thereby increasing the 
space between the keys. 

The keys are divided into 
groups according to func- 
tions. The groups are sepa- 
rated by size, value contrast, color, and placement 
of nomenclature. The numeric keys, which are most 
frequently used, are larger and have the strongest 
value contrast. They have their nomenclature di- 
rectly on the keys. The next group of keys accord- 
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ing to frequency of use are identified by their blue 
color. The ENTERt key and arithmetic keys are 
separated within this group by placement of nomen- 
clature on the keys. The less frequently used keys 
have the least value contrast, and the nomenclature 
is placed on the panel just above the key. 

The keys have an over-center or snap feel when 
they are pressed. This comes from special spring 
contact developed by HP. The electrical and me- 
chanical parts of the keyboard are less than "h inch 
high. 

The spring contacts are mounted on a printed- 
circuit board along with several other components 
including the fifteen-digit I.ED display. The display 
is tilted toward the operator for optimum viewing. 

A second printed-circuit board carries the ma- 
jority of the electronic components. The two boards 
are attached by a series of pin connectors. 

The external package was developed from a hu- 
man-engineering approach, with aesthetic appeal of 
major importance. The sculptural wedge shape per- 
mits the calculator to be comfortably held in the 
palm of one hand. It also allows the product to slide 
easily into a pocket. The keyboard and display slope 
upwards for a better viewing angle in desk-top use. 
The sculptural sides visually break up the total 
mass of the package. The top half of the case is 
highlighted while the bottom half is in shadow. 
This gives the product the appearance that it is 
thinner than it actually is. The product appears to 
be floating when viewed from a normal operating 
position in desk-top use. The use of textures that 
complement each other contribute significantly to 
the overall finesse and appearance. The texture on 
the case provides a non-slip surface, important 
when the calculator is being hand-held. 

Viewed from the bottom, the calculator retains 
a clean appearance. There are no exposed screws 
or other fastening devices. This is aesthetically im- 
portant to a product that is hand-held and viewed 
from all sides. 

The polyvinyl-chloride feet prevent the calcula- 
tor from sliding during desk-top use. The rear feet 
also serve as the battery-door latches, aiding in the 
overall cleanness of the product. 

The HP-35 couldn't have been developed without 
an outstanding working relationship between labo- 
ratory, industrial design, manufacturing, and tool- 
ing. Everyone involved in the project shared a 
common desire to retain the original size and shape, 
and many innovative engineering concepts resulted. 
Many of the problems encountered during develop- 
ment could have been easily solved by using more 
conventional methods, but the result would have 



been a larger package. Now that the calculator is a 
reality, everyone feels that the extra efforts re- 
quired were worthwhile. J 
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SPECIFICATIONS 

HP-35 
Pocket Calculator 

FUNCTIONS: 

ARITHMETIC: Add, Subtract. Multiply. Divide and Square Root. 
TRIGONOMETRIC: Sin », Cos x. Tan k, Are Sin ». Arc Cos x, Arc 
Tan x. 

LOGARITHMIC: Log | >. Log. >. and 8' 

OTHER FUNCTIONS: X', 1/x, r and data jlorage and poallloning 
keys 

SPEED OF OPERATION (typical): 

Add. Subtract 60 mlllleecondi 

Multiply. Divide 100 milliseconds 

Square Root 110 milliseconds 

Logarithmic & Exponential 200 milliseconds 

x> 400 milliseconds 

Trigonometric 500 milliseconds 

POWER: 

ac: 115 or 230 V i:10%, SO to 60 Hi. S watts. 

Battery: 500 mW derived from Nickel-Cadmium rechargeablo 
Battery Pack. 
WEIGHT: 

Calculator: 9 oz,, Recherger 5 oz 

Shipping weight: approx. 2 lbs. 
TEMPERATURE RANGE: 

OPERATING OX to 40X (32"F to 104 - F) 

STORAGE: -40X to + 55'C |-40*F to 131'F). 
PRICE IN USA: $395, including battery pack. 115/230 V adapter/ 

recharger carrying case travel case, name tags, operating manual 
Note: At press time, orders for the HP-35 have exceeded expecta- 
tions to such an extent that a waiting list has been established 
Deliveries should Improve In the next few weeks. 

MANUFACTURING DIVISION: ADVANCED PRODUCTS DEPARTMENT 

10900 Wolfe Road 
Cupertino, California 95014 
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New Capabilities in Digital 
Low-Frequency 
Spectrum Analysis 

A new Fourier analyzer and two fast-transform periph- 
erals adapt to a wide range of applications 

By Stephan G. Cline and Norman D. Marschke 



SPECTRUM ANALYSIS, a powerful technique 
of signal and system analysis, is the process of 
resolving a physical waveform into sinusoidal com- 
ponents having appropriate phase relationships. 
That this can be done was shown by J.B.J. Fourier 
in the 18th century. The reason for its great utility 
is that a simple statement of spectrum content or 
bandwidth often makes possible engineering judg- 
ments about a problem that couldn't be made on the 
basis of observations of the waveform itself. 

In spite of its usefulness, applications of spec- 
trum analysis were concentrated for many years in 
the radio and electronic sciences, where the fre- 
quencies encountered are in the megahertz range 
and higher. Now, however, the use of spectrum 
analysis at low and moderate frequencies is on the 
increase, partly because many more people realize 
its power as an engineering tool, but mainly because 
accurate and stable instrumentation has recently 
become available for this frequency range. Prob- 
lems in mechanical vibrations, underwater sound, 
biomedicine, communications, machine mainte- 
nance, feedback control systems, speech, and acous- 
tics are now yielding to spectrum analysis. The in- 
struments being used are mostly digital analyzers 
which compute spectra by means of one form or 
another of the fast-Fourier-transform (FFT) algo- 
rithm. 1 

Within the last year, faster computers and hard- 
ware FFT processors have significantly increased 
the speed of Fourier analysis, to the point where 
real-time spectrum analysis to beyond 10 kHz can 
be done using entirely digital instruments. Among 
Hewlett-Packard's second-generation Fourier-anal- 
ysis instruments are a new Fourier Analyzer, a new 
plug-in FFT Arithmetic Unit for HP computers, and 
a new peripheral Fast Fourier Processor. 



Why Digital? 

Why did the use of spectrum analysis at low fre- 
quencies have to wait for the development of digital 
instruments before becoming widespread? The an- 
swer can be seen by comparing digital analyzers 
wilh previously available instruments. 

The simplest way to resolve a waveform into its 
spectral components is to use a narrow bandpass 
filter and a voltmeter. Measuring the portion of the 
signal passed by the filter gives an approximation 
to the size of the sinusoidal component at the filter 
center frequency. If the filter is tunable, spectrum 
analysis can be carried out over a range of frequen- 
cies. This is the principle of the wave analyzer, the 
most elementary spectrum analyzer. 

The wave analyzer has the basic limitation of all 
spectrum analyzers. The bandwidth of the filter, 
and therefore its resolution, determines its response 
time. A period of time approximately equal to the 
reciprocal of the analyzer passband must be spent at 
each frequency to get an accurate estimate of the 
spectral content at that frequency. This isn't a seri- 
ous problem at frequencies in the megahertz range 
and higher, but at low frequencies, where the reso- 
lution needed may be on the order of one hertz, a 
single-frequency measurement may take as much as 
one second. 

If the signal being analyzed is random, each read- 
ing must be averaged over an even longer period of 
time to get a statistically stable or smooth result. A 
measurement made in this way may correctly meas- 
ure the spectrum at 1.0 kHz while missing a short 
term event at 1.5 kHz. 

For these reasons, wave analysis is suitable only 
when the signal being analyzed is stationary and 
when comparatively long periods of time are avail- 
able to sweep over the band of interest. Heterodyne 
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Fig. f a Power spectra ol random noise illustrate resolution, dynamic range, and flexibility 
ot digital analyzers, fa) Spectrum computed trom one 102 4 ms record (b) Average ol 200 
records has same 512-channel resolution as single record. Dynamic range is BO dB. Note 
1.2-kHz component at about —70 dB which is not visible in unaveraged spectrum. 



techniques which translate the measurement to 
higher frequencies make for better accuracy, better 
display, and less cost, but they don't solve the time- 
bandwidth limitation at low frequencies. Spectrum 
analyzers of this type are still basically single-filter 
wave analyzers. 

The measurement time-bandwidth problem can 
only be alleviated if some form of parallel filtering 
is employed. The most direct way to build a spec- 
trum analyzer that is real-time at all frequencies is 
to use a set of parallel filters. Such a method is most 
effectiv e when filters of constant percentage band- 
width, such as Va octave or full octave, are em- 
ployed. Octave-band analyzers of this type are very 
useful and effective instruments. 

Speed-Up Analyzer 

However, when higher resolution or constant- 
bandwidth filters are needed, this approach has 
proved impractical. The nearest solution to parallel 
filtering using constant-bandwidth, narrow-band 
filters is a digital-analog hybrid approach called the 
speed-up analyzer. In this analyzer an analog-to- 
digital converter samples the waveform to be ana- 
lyzed at a rate F s and loads a memory with a record 
of N samples over a period of time T = N/Fs. 
The data gathered in this way is analyzed by re- 
circulating it in the memory at a much higher rate. 
The data time scale is thus reduced and the speeded 
up digital data is applied to a digital-to-analog con- 
verter, resulting in a new signal with scaled-up fre- 
quency, This new signal is then analyzed by the 



equivalent of a high-frequency wave analyzer. At 
the higher frequency the analyzer sweep time is 
greatly reduced and near-real-time analysis can be 
carried out*. It isn't quite real-time, though, because 
there's a 20-to-100-millisecond gap between the 
processing of one T-second record and the process- 
ing of the next, so some of the signal information is 
lost. 

All-Digital Analyzer 

An important limitation on the performance of 
speed-up analyzers is the loss of dynamic range and 
resolution introduced by the analog circuitry of the 
final wave analyzer. Better dynamic range and reso- 
lution and greater flexibility in data processing can 
be obtained with a purely digital analyzer. 

Like the speed-up analyzer, the all-digital ana- 
lyzer samples the input signal at a rate F s for T sec- 
onds and loads N = TFs sample values into its 
memory. Certain measurement limits imposed by 
the sampling process are shared by the digital 
analyzer and the speedup analyzer. First, the highest 
frequency that can be analyzed without ambiguity 
is 1/2 F s . This maximum frequency is the Nyquist 
folding frequency; above it. aliasing error occurs-. 
Second, the resolution is determined by the record 
length. A record length of T seconds implies that 
the frequency resolution is Af = l/T. 

This smallest frequency element is a result of the 
fact that for a record of length T seconds no discrete 

■Real-time analysis it Defined ne'e as ft analysis of contiguous records of aata *!tn 
no gsos. 
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Fourier transform can be defined with a resolution 
narrower than Af Hz. It has nothing to do with any 
other hardware constraint such as filter band- 
widths. 

The digital analyzer transforms the N digitized 
values stored in its input memory to a set of N/2 
Fourier coefficients spaced 1/T Hz apart. The re- 
sult of this Fourier transform is a set of real and 
imaginary coefficients, which may be used directly 
or further processed to obtain the magnitude and 
phase of the spectrum, or the power (magnitude 
squared) spectrum, or many other functions. 

The measurement time-bandwidth limitation of 
wave analyzers isn't a problem with digital ana- 
lyzers. From a single T-second sample record of N 
samples, the digital analyzer computes a complete 
spectrum equivalent to N/2 wave-analyzer measure- 
ments each 1 T Hz apart. The same spectrum would 
take at least N/2 times as long to measure with a 
wave analyzer. 

If the data being analyzed is random or has a ran- 
dom component, so that several separate sample 
records must be averaged to smooth the result, this 
is easily done in digital spectrum analysis by sum- 
ming the results of one sample record into the ac- 
cumulated result from past samples. One of the 
advantages of digital computation of the spectrum 
is that averaging may be done digitally without the 
use of analog integrators, which have limited dy- 
namic range and linearity. 

An Example 

Fig. 1 illustrates some of the capabilities of the 
latest digital analyzers. It shows the power spec- 
trum of a band-limited random noise signal com- 
puted in real time. Each record used to compute 
the spectrum is 102.4 milliseconds long and consists 
of 1024 samples of the noise signal. 

Fig. 1(a) is the spectrum computed from one rec- 
ord, and the variance caused by the randomness of 
the signal is evident. Fig. 1(b) is the average of the 
spectra computed from 200 sample records. Besides 
being much smoother, this spectrum reveals the 
presence of a periodic component at 1.2 kHz. These 
results, computed from only 20.48 seconds of data, 
are good examples of the resolution, dynamic range, 
and processing flexibility of a digital analyzer. 

New Fourier Analyzer 

Fig. 1 was produced by the new Model 5451A 
Fourier Analyzer, Fig. 2. Like earlier HP Fourier 
analyzers^, the new analyzer is a keyboard-con- 
trolled computer-based system capable of sampling 
one or more input signals and computing power 




Fig. 2. New Model 54S1A Fourier Analyzer, a minicom- 
puter-based, keyboard controlled, dual-channel analyzer, 
can do spectrum analysis, correlation, coherence analysis 
and many other complex data manipulations. New capabil- 
ities are 80-dB dynamic range, relocatable software, but- 
tered input, and a faster computer. 

spectra, cross power spectra, transfer functions, 
coherence functions, correlation functions, and 
many other functions. These computations and 
other complex data manipulations are performed in 
the computer by software routines which are ini- 
tiated by pushing buttons on the keyboard. The user 
can choose to add his own routines to the software 
supplied. He can also use the integral minicomputer 
separately as a general purpose computer. Fig. 3 is 
a block diagram of the Fourier analyzer system. 

Differences between the new 5451A and earlier 
models are a new faster computer, Model 2100A*. 
greater dynamic range (80 dB for both amplitude 
and power spectrum), relocatable software which 
simplifies adding and modifying user-written pro- 
grams, and a buffered input mode which improves 
processing efficiency. 

Fig. 4 shows the difference between the new buf- 
fered input mode and the ordinary input mode. Or- 
dinarily, an input record is read into the computer 
memory and then processed, and while the power 
spectrum is being computed no new input data is 
collected. Without its optional hardware proces- 
sors, the 5451A Fourier Analyzer may take as long 
as 1.5 seconds to complete these processing steps 
for a 1024-point record. Considerable improvement 
is realized when the buffered input mode of opera- 
tion is used. In this mode data is read into a buffer 
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Fig. 3. Analog-to-digital converter samples input signals 
al rales lo 50 kHz Keyboard controls data manipulations 
done by software in basic 5451 A Fourier Analyzer. Fast- 
processor options add hardware transform capability. 

block while the previous record is being processed. 
With this parallel input-while-processing scheme it 
is possible to achieve real-time spectrum analysis 
to .'lf)0 Hz using a 5451 A Fourier Analyzer with only 
a software Fourier transform. 

Faster Fourier Analysis 

Further improvements in processing speed are 
realized with special-purpose hardware which 
speeds up the required computations. Either the 



T)470A Fast Fourier Processor or the 5471 A Fast 
Fourier Transform Arithmetic Unit can be added 
as options to the 5451 A Fourier Analyzer for this 
purpose. 

Increases in processing speed increase the ana- 
lyzer's upper frequency limit for real-time analysis. 
Increased speed also saves machine time, which is 
important in test situations of limited duration, such 
as vibration testing of a spacecraft which might be 
fatigued in a long test. 

The 5470A Fast Fourier Processor (Fig. 5) is an 
external peripheral which has its own data memory 
and function control memory. It takes only 15 milli- 
seconds to compute a Fourier transform for a 1024- 
point block size, compared to one second for the 
standard 5451 A Fourier Analyzer. 

The four-card 5471A FFT Arithmetic Unit (Fig. 6) 
fits into four contiguous I/O slots in the computer 
and speeds up computations while continuing to 
use the computer for data storage and function con- 
trol. With the 5471 A a 1024-point transform can be 
computed in 100 ms. 

The software drivers for the 5470A and 5471 A 
are written so analyzer operation is exactly the 
same as the standard 5451 A Fourier Analyzer, ex- 
cept that the execution times are much faster. 

The following array operations are performed by 
the 5470A or the 5471A options in the 5451A. 

• Forward and inverse fast Fourier transform. 

• Time-domain weighting by Hanning. 

• Self-conjugate block multiplication in single or 
double precision. 

• Block addition/subtraction and scaling in 
single or double precision. 



Unbuffered Input Mode: 

lo^ita ^tak' I Process Record 1 I ln P ul Record 2 j Process Record 2 Input Record 3 



Time 



Buffered Input Mode: 



'""to Buffer 11 Input Record 2 | Input Record 3 Buffer Operations 

Process Record 1 1 | Process Record 2 1 | Processor Operations 



Transfer Data 
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to Data Block- 



Fourier 
Transform 



Self Conjugate 
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Data Block Locations: 

5451A: In computer memory. 
5451A/5470A: In 5470A 
5451A/5471A: In computer memory. 



Fig. 4. Timing diagrams lor 
computation of power spectrum 
average. New buttered input 
mode is capability ol 5451 A 
Fourier Analyzer with or without 
lasi-orocessor ootions. Parallel 
input-while-processing ooera- 
tion in buttered mode increases 
system etliciency. 
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Multiple array operations are used to perform com- 
plete spectrum analysis functions such as: 

• Auto and cross power spectrum average. 

• Auto and cross correlation. 

• Convolution. 

• Transfer and coherence functions with the aid 
of division routines in the host computer. 



Fast Fourier Processor 

The 5470A Fast Fourier Processor is a digital 
processor designed to perform the Fourier analysis 
routines rapidly and at reasonable cost. As a pe- 
ripheral to the 5451 A system, the 5470A operates 
under functional control of the host computer which 
provides it data, microprograms, and execution 
command parameters. 




Fig. S. New Model 5470 A Fast Fourier Processor used with 
5451 A Fourier Analyzer lowers time tor 1024-point trans- 
form to 15 milliseconds Processor can also add Fourier - 
transform capability to other computer systems. 



Nearly all data-array manipulations normally 
needed for spectrum analysis are performed in the 
Fast Fourier Processor. It has either a 4096-word 
data memory or an optional 8192-word data mem- 
ory, and data blocks may range in size from 64 to 
2048 words (4096 optional). 

When a power spectrum average is being com- 
puted, data flow in the combined 5451 A/5470A sys- 
tem is similar to that in the analyzer alone. After 
the buffer in the computer is filled with input data 
from the analog-to-digital converter, the data is 
transferred to the data memory in the 5470A. For 
each data record a time-to-frequency Fourier trans- 
form, a self-conjugate multiply, and an addition into 
the accumulation block are successively performed 
in the 5470A. While these operations are in progress 
the computer is free to load the next data record 
into the buffer section of its memory from the ana- 
log-to-digital converter, and send output to the dis- 
play if time permits. 

All block operations on data in the 5470A are 
controlled by microprograms from the computer. 
For each operation a microprogram from the soft- 
ware driver is loaded into the processor's function 
control memory and then executed. Each user com- 
mand may call for multiple 5470A operations. 

The Fast Fourier Processor is linked to the com- 
puter in the 5451A Fourier Analyzer by two I/O 
channels which use standard 16-bit TTL micro- 
circuit interface cards. Since the interface cards are 
standard, the fast processor doesn't have to be used 
with a Fourier analyzer: it can add Fourier-analysis 
capability to other computer systems as well. 
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The Where and How 
of Fourier Analyzers 

Here are examples of how HP Fourier Analyzers are being 
used in various locations around the world: 
Massachusetts. Used aboard a research ship for sea floor 
exploration in search of oil fields. 

Canada Underwater Acoustics Extensive use of coherence 
and transfer function analysis to determine optimum sonar 
bubble shape. 

Germany Acoustic and vibration analysis of linear servo 
systems. 

California Civil Engineering. Used to study dynamic charac- 
teristics of large structures, such as resonant frequencies 
and damping factors, to determine how they will with- 
stand earthquakes. 

Sweden. Biomedical. Electromyography. Used in studying 
types of signals transmitted on various types of nerves for 
designing automatic prostheses 

Michigan. Vibration analysis on vehicle drive trains in pro- 
duction environment to determine correct/incorrect opera- 
tion of gears and other parts. 

Iowa. Vibration analysis on tractors to aid In design of drive 
train and frame components. 

Belgium. Vibration analysis on machinery used in the ap- 
plication of emulsion to film to determine effect of vibrations 
on emulsion. 

New York. Vibration analysis on cutting tools to determine 
wear characteristics. 

France. Vibration studies of helicopters aimed at noise re- 
duction. 

Netherlands. Civil Engineering Analysis of vibrations of 
tunnels, dikes and other large structures to determine 
dynamic characteristics. 

Michigan. Vibration analysis of production power tools to 
determine wear and efficiency. 

Michigan Vibration analysis of axles and axle housings to 
improve design and for production testing. 
Michigan. Vehicle crash studies Used to design automobile 
components such as steering wheels. 
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Fig. 6. New Model 5471A Fast Fourier Transform Arith- 
metic Unit plugs into tour I/O slots in HP computers. 1024- 
point transform time is 160 milliseconds. 



Fast Fourier Transform Arithmetic Unit 

The 5471 A Fast Fourier Transform Arithmetic 
Unit speeds up by a factor of 6 or more many of the 
operations required in Fourier analysis. The FFT 
Arithmetic Unit is a hardware-software package. 
The hardware consists of four printed-circuit 
boards with plug into any four contiguous I/O slots 
in a Hewlett-Packard 2100-series computer. 

The software consists of a driver written in HP 
assembly language. The driver has been incorpo- 
rated into the Fourier analyzer software so to the 
user the operation is the same as the standard ana- 
lyzer except for the faster speed. Another assembly 
language driver allows the 5471A to be used with a 
stand-alone computer by those users who need 
Fourier arithmetic capability but don't want a com- 
plete Fourier analyzer system. 

Data is stored in the computer, not in the arith- 
metic-unit hardware. The 5471A memory is simply 
a small scratchpad. Memory addressing and read/ 
writing are done by software while the arithmetic 
operations are done by the 5471A hardware. 



SPECIFICATIONS 
HP Model 5451 A 
Fourier Analyzer 

ANALOG INPUT 

AMPLITUDE RANGE: 0.1 V to 10 V lull scale. 
RESOLUTION: 10 oils with sign. 

DISPLAY UNIT 

Data may be displayed on the 8 x 10 cm CRT display or output 
to a plotter or remote display: 

AMPLITUDE SCALE: Data In memory Is automatically scaled to give 
a maximum on-screen calibrated display. 
LINEAR DISPLAY RANGE: ~* divisions with scale (actor ranging 

from 1 X 10-:" to S x 10-«> In step* ol 1. 2, 5, and 10. 
LOG DISPLAY RANGE: 4 decades with a scale factor ranging 
from 0 to 3:998 dB. 

DIGITAL ACCURACY AND RESOLUTION 

All calculations use floating point arithmetic on a block basis. 
Data overflow does not occur. Amplitude resolution Is 1 part In 
16,000 worst case. 

DATA WORD SIZE: 16-blt real and 18-blt Imaginary or 16-blt mag- 
nitude and 16-bit phase. 32 bits real Is preserved for auto power 
spectra. Division, addition, and subtraction operations performed 
In 16 blta or 32 bits depending on data 

TRANSFORM ACCURACY: 0 1% of full scale worst case error 
during forward or Inverse calculation. 

TRANSFORM SIZE: From 64 to 4096 time points. 

TRANSFORM SPEED: 1 s for 1024 time points. 

SPECTRAL RESOLUTION 

The element of spectrsl resolution la the frequency channel width 
— the maximum frequency divided by Vj the data olock size 

MAXIMUM FREQUENCY: 25 kHz single channel; 10 kHz dual chan- 
nel Adjustable In steps of 1. 2.5. 5, and 10 down to 0.1 Hz. 

DYNAMIC RANGE: 80 dB for voltage magnitude spectrum end auto 
power spectrum. 

POWER SOURCE: 115/230 volta rr10%, 50/60 Hz. 1400 watts typical 
for baalc system (excluding teleprinter). 

ENVIRONMENTAL CONDITIONS: O'C to 55*C. 

Relative Humidity: To 95% at 40'C (104'FI. 
PRICE IN USA: 5451A System prices star; at S35.000 



HP Model 5470A 
Fast Fourier Processor 

TRANSFORM BLOCK SIZE: 64 :o 2048 real points standard. Up to 

4096 points optional. 
DATA MEMORY: 

4096 16-blt words, optional 8K available. 
SINGLE PRECISION OPERATIONS: 

Fest Fourier transforms — forward and Inverse. 

Interval centered Hannlng. 

Self multiply: real, complex, and con|ugete. 

Cross multiply: real, complex, and conjugate. 

Scale and add. 
0OU8LE PRECISION OPERATIONS: 

Self conjugate multiply. 

Scale and add. 

FOURIER TRANSFORM TIME (based on a 1024 real point block 

size): 15.2 ma. 
ENVIRONMENTAL: O'C to 55'C. 
PRICE IN USA: 5470A p'ices start at S25.000. 

HP Model 5471A 
Fast Fourier Transform Arithmetic Unit 

ARRAY SIZES 

TRANSFORM BLOCK SIZE: 64 to 2048 points In powers of 2. 
ADDITION BLOCK SIZE: 64 to 2048 points In powers of 2. 
SINGLE PRECISION OPERATIONS 
Forward and inverse fast Fourier transforms. 
Interval centered Hannlng. 

Real, complex, and complex conjugate multiply. 

Block scaling and block addition. 
DOUBLE PRECISION OPERATIONS 

Self complex conjugate multiply. 

Block scaling and block addition, 
FOURIER TRANSFORM TIME (based on 1024 word block size): 

160 mi with 2100A Computer. 
ENVIRONMENTAL: 0' to 55'C. 
PRICE IN USA: 5471A prices start at $4500. 
MANUFACTURING DIVISION: SANTA CLARA DIVISION 

5301 Stevens Creek Boulevard 
Santa Clara, California 95050 
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The 5471A must be microprogrammed for each 
type of arithmetic operation (i.e., Fourier transform, 
complex multiply, double-precision self-conjugate- 
complex multiply, Hanning, etc.). For each opera- 
tion a microprogram from the software driver is 
loaded into the program memory on the 5471 A pro- 
gram board and then the 5471A is ready to start. 

The design philosophy of the 5471A is analogous 
to the concept of impedance matching for optimal 
power transfer. The 5471 A matches times instead 
of impedance. It makes the time required for the 
hardware to perform its arithmetic tasks approxi- 
mately equal to the time required for the software 
to output 16-bit data words to the 5471 A, generate 
the next data addresses, and then input IB-bit data 
words from the 5471A to the computer. The arith- 
metic operations consume approximately 85% of 
the machine time required to do a Fourier transform. 
The 5471A hardware actually processes the data 
fasler than the computer can calculate the next 
addresses. Thus the 5471 A is able to speed up the 
transform time by a factor of 6 to 7. 
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